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ABSTRACT 

We investigate the claim that the ratio /? of radiation pressure force to gravitational force on a 
dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain 
porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we 
find that the maximum value of f3 is almost independent of grain porosity, but for small (< 0.3/tm) 
grains, (3 actually decreases with increasing porosity. We also investigate the effect of metallic iron 
and amorphous carbon inclusions in the dust grains and find that while these inclusions do increase 
the radiation pressure cross-section, (3 remains below unity for grains with 3 pg of silicate material. 
These results affect the interpretation of the grain trajectories estimated from the Stardust mission, 
which were modeled assuming j3 values exceeding one. We find that radiation pressure effects are not 
large enough for particles Orion and Hylabrook captured by Stardust to be of interstellar origin given 
their reported impact velocities. We also consider the effects of solar radiation on transverse velocities 
and grain spin, and show that radiation pressure introduces both transverse velocities and equatorial 
spin velocities of several hundred meters per second for incoming interstellar grains at 2 AU. These 
transverse velocities are not important for modeling trajectories, but such spin rates may result in 
centrifugal disruption of aggregates. 


1. INTRODUCTION 

A longstanding question in the field of solar system 
dust grain dynamics has been the value of the ratio (3 of 
radiation pressure force to gravitational force on a dust 
grain. Because both forces scale with the inverse square 
of the distance to the Sun, this ratio is independent of 
location in the solar system. There has been widespread 
agreement for several decades that as a function of grain 
size, f3 peaks somewhere around one for grains with radii 
of a few tenths of a /am, but it has been a matter of de¬ 
bate on which side of unity the peak lies. This is of some 
practical significance, since a value of (3 greater than one 
for a certain grain size means not only the absence of or¬ 
biting grains of that size, but also a paraboloidal region 
of space wh ere n o i nterstellar grains of that size can pen¬ 
etrate (Landgraf et al. 1999b). The Ulysses and Galileo 
missions found a deficit of grains in the mass range from 
0.01 to 0.3 pg (corresponding to radii of a few tenths of a 
/mi) inside of 4 AU in our solar system, which was inter¬ 
preted as being due to radiation pressure excl uding t hose 
grains f rom the inner solar system (e.g., Landgraf et al. 
1999a). Kimura et al. (2003) found that the Ulysses data 
a re best tit assuming p < 1 for grains larger than 0.3 pg. 


Burns et al. (1979) used Mie theory to calculate f3 for 
spherical grams composed of a variety of materials. They 
found /3 values peaking around 0.6 for quartz, but well 
above unity for certain other materials such as graphite, 
iron and ma gn e tite. T here has been speculation (e.g. 


Landgraf et al. (1999a}) that porous grains may have a 
significantly higher value of (3. This is intuitively ap¬ 
pealing, sinc e the geometric c ross-section per unit mass 
is increased. Saija et al. (2003) studied the effect of grain 
structure on small (tens of nanometers) grains, and found 
less compact grains to have lower f3 values. We extend 
these computations to grain sizes as large as a e ff = 1/un, 


where a e g = (3U/(47 t)) 1 ^ 3 is the volume-equivalent ra- 
d ius for a grain wi th sol id volume V. 


Kimura et al. (2002) studied radiation pressure on 
grams formed through ballistic aggregation. Their cal¬ 
culations were limited to grains with volume-equivalent 
radii a e ff < 0.20/mr. They used a dielectric function 
appropriate for obsidian, which has very little absorp¬ 
tion. In their calculations they found that their aggregate 
structures have /3 values nearly an order of magnitude 
lower than those of a solid sphere for a e ff = 0.1/im, with (3 
dropping to very low values as a e R shrinks below 0.1 /mi. 
In this paper, using a dielectric function appropriate for 
astrosilicates, we find more modest differences between 
our aggregate models and the single-sphere case, and we 
find that f3 goes to a constant value for a e g < 0.1/rm, as 
expected when A a e g and absorption dominates over 
scattering. 

The results of the Stardust mission provide new reason 
to pin down values of /3 for large r grains with com plicated 
geometries. The Stardust team (Westphal et al. 2014) has 
identified three a e g « 0.6/rni captured particles which 
they believe to be of interstellar origin. The low impact 
speeds necessary for particles to survive the collection 
process require /? > 1 so that solar radiation can deceler- 
ate the incoming particles prior to impact (Sterken et al. 
This conflicts with the results from Mie theory 
r show /3 peaking around 0.8 for silicate grains. 


20141. 

whic 


In this paper, we use accurate scattering calculations 
for model dust grains to examine the claim that more 
porous dust grains have values of [3 well above th ose fo r 
solid spheres. An alysis of the Stardust particles (West¬ 


phal et al. 2014) shows complicated compositions, in¬ 
cluding some containing “Fe-bearing phases” said to be 
consistent with metallic iron. To address this, we also 
calculate f3 for grains with metallic iron inclusions. 
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TABLE 1 
Filling factors 


Fig. 1.— Representatives of each class of grain for N = 32. The 
vectors a; represent the principal axes of the grain. 

2. TARGETS 

The grain parameters that determine /3 are grain size, 
shape, material density, and dielectric function. We char¬ 
acterize size with an effective radius a e g. This is the 
radius of a solid sphere with equal mass. For grain 
sh apes, we us e reali zations of the random aggregates 
of|Shen et ah ([2008]), taken from http://www. astro. 
princeton.edu/~draine/agglom.html These model 
grains are formed through ballistic aggregation, and Shen 
et al. provide three algorithms which yield grains with 
different porosities. In their most porous model, de¬ 
noted “BA,” a grain starts as a single sphere. Addi¬ 
tional spheres approach on randomly oriented trajecto¬ 
ries, and if they impact the grain, they stick to the first 
point of contact. In their intermediate density aggregates 
(denoted “BAM1”), once each new sphere has impacted 
the target, it rolls along the surface to contact the next 
nearest sphere that is part of the grain. In their most 
compact structures (“BAM2”), this happens once again 
so each new sphere after the third ends up touching 3 
spheres already part of the target before approach of the 
next impactor. Shen et al. provide realizations of these 
structures for numbers of spheres between N = 2 3 and 
2 16 . In Figure[l]we show visualizations of one grain from 
each class for the realizations with 32 spheres. 

For most of these calculations we use the realizations 
with N = 32, although we verify that results are not 
dramatically changed by considering realizations with 
larger N. Shen et al. provide average filling factors / 
and porosities 1 — / for these structures. They deter¬ 
mine / by finding the uniform density ellipsoid with the 
same mass and moment of inertia tensor as the grain. 
The filling factor / is then the volume of the constituent 
spheres divided by the volume of this equivalent ellipsoid. 

The filling factor / varies substantially between the dif¬ 
ferent grain models considered in this paper (from 0.15 
for the BA, N = 256 model to 0.55 for the BAM2, N 
= 32 model), but is fairly uniform between different re¬ 
alizations of models, all models having standard devia- 
tion o f filling factor ctf on the order of 10% (Shen et al. 
2008). Mean fillingfactors and standard deviations are 
presented in Table |T) 

For our grain composition, we use the die lectric 
functi on estimated for amorphous MgFeSi04 in Draine| 
(2003). This silicate material is assumed to have a den- 
sity p s = 3.8 g cm -3 , intermediate between forsterite 
Mg 2 Si04 (p s = 3.21 g cm -3 ) and fayalite Fe 2 SiC>4 (p s 
= 4.39 g cm -3 ). We also consider grains with inclusions 
of metallic iron and amorphous carbon. The iron inclu- 
si ons have a dielectri c funct ion described in appendix B 
of |Draine fe Hensley (2013) and density appropriate for 
metallic iron (7.87 g cm -3 ). We consider 3 carbon dielec¬ 


Grain Class 

N 

/ 

°7 

°7 /(f) 

BAM2 

32 

0.55 

0.031 

0.056 

BAM2 

256 

0.42 

0.026 

0.061 

BAM1 

32 

0.36 

0.034 

0.094 

BAM1 

256 

0.26 

0.018 

0.068 

BA 

32 

0.20 

0.031 

0.155 

BA 

256 

0.15 

0.013 

0.086 


tric functions taken from Zubko et al. ( 1996|), and assume 
a carbon density of 2.0 g cm -3 (Robertson] 1986). 


We must also specify the orientation of the grain with 
respect to the direction of incident radiation (i.e. the 
Sun). This is a complicated subject, requiring averag¬ 
ing over a set of orientations depending on the spin state 
of the grain. Fortunately, we find that orientation only 
affects /3 at the 10% level. This is discussed further in 
Section 15.21 

3. RADIATION PRESSURE CODES 

In this paper we use two different codes to calculate 
radiation pressure on dust grains. First, we use the Mul¬ 
tiple Sphere T Matrix (MSTM ) code written by Daniel 

This code calculates ab- 


Mackowski (Mackowski 2013). 
sorption and scattering by a collection of spheres. 


The 


algorithm is based on decomposition of the scattered 


around each of the sc 

iheres 

(Mackowski 1994 

Mackowski 

& Mishchenko 

1996 

2011 

. The MSTM method is in 


ily truncated at a finite number of spherical harmonics 
for each constituent sphere (just as in Mie theory calcu¬ 
lations for a single sphere). 

We also used the d i screte dipole code DDSCAT 
(Draine & Flatau 1994 20131. It gives, among other 


things, the components of the mean scattered radiation 
field in the transverse directions. For the radial radiation 
pressure, where we can do a direct comparison between 
the DDSCAT code and the MSTM code, we found agree¬ 
ment to within a few percent. The MSTM code is not 
configured to easily give the transverse components of 
the scattered field, or to calculate torques on grains, so 
for this we used DDSCAT, even though it is substantially 
slower. 

4. CALCULATION OF /3 

For a given grain geometry, incident wavelength A, and 
dielectric function, we calculate Q a b s , Q sea and (cosd sca ), 
where 9 sca is the angle a scattered photon makes with the 
direction of the incident beam, <5 a b s = Cabs/OTCQj), and 
Qsca = C sca /(Tral s ). C^s and C sca are the cross-sections 
for absorption and scattering respectively. The radial 
radiation pressure force due to sunlight is given by 


Fr a d — 


7ra^ ff Q pr (A)F 1 (A)(iA 


where 


Qpr — Qabs Qsca(l (cOS$ sca )), 


(i) 


( 2 ) 


F(X)dX is the flux of sunlight in [A, A + dX\. and c is the 
speed of light. We approximate the Sun by a blackbody 
with temperature T 0 = 5777 K ( CoxpO OO). 
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Fig. 2.— Calculations of (Q a bs)©, (Qsca)© and ((cos0 sca )) Q 
from both MSTM and DDSCAT. This is for a particular orienta¬ 
tion of a particular instance of a grain of the BAM1 class with 32 
spheres (BAM1.32.9). The MSTM and DDSCAT results are nearly 
indistinguishable. 


We define (x)@ as the quantity x averaged over the 
solar spectrum. In the case of the radiation pressure 
efficiency Q pr , 


, f 0 °° Q pr F(X)dX 

Wprio- joc F(A)dA 


Q pr (X)B x dX , 


(3) 

where B\ is the Planck function corresponding to T 0 , 
and a is the Stefan-Boltzmann constant. 

In Figure [2] we show a plot of (Q abs } 0 , (Qsca)o and 
((cos$s C a)}© as a function of a e ff.n ((cos 9 scs ))q is the 
mean value of (cos 9 sca ) weighted By F(X)Q sca . We plot 
the results from both the MSTM code and DDSCAT. 
We see that the scattering cross-section is dominant over 
absorption except for very small sizes. Their relative 
effects on radiation pressure are more equal because of 
the preference for strong forward-scattering when Q sca > 
1. At distance R from the Sun, 


F( X) = ^ B x , (4) 

where i? 0 is the radius of the Sun. Using equations ([I]), 
([3]) and Q, we find 

p (n y ™lsP T Q R l, n \ 

Trad(aeffj — -- {Q pr/0- (5) 


2 See http://www.astro.princeton.edu/~draine/SD2016.html 
for additional details on computation of the results shown in Figure 
2 . 


Fig. 3.— The ratio Q P r(A)/a e ff for different values of a e ff for a 
particular orientation of a particular TV = 32 grain of the BAM1 
class (BAM1.32.1). Also plotted is the normalized solar spectrum 
per unit ln(A) (black curve). The black dots are the A i at which 
we sample the solar spectrum. The product of these two curves is 
proportional to the grain acceleration per unit interval in ln(A). 

Dividing by the gravitational force on the grain 




Fr&d 

F 

x grav 


3 Lq (Qpr)@ 

167tcGMq p s a e g ’ 


( 6 ) 


where M 0 is the mass of the Sun and G is the gravi¬ 
tational constant. To calculate (Q pr )o numerically, we 
approximate F(X) as a sum of N\ delta functions, each 
with -p- of the power of the Sun, spaced such that the 
i th delta function is located at A i such that 


F(X)dX = <71% R& r 


i- 1/2 


and 


R 2 L N x 
1 N > 


i = l,...,N x 


(Qpr)© ~ jy Q P r(Xj). 


»= 1 


(7) 

( 8 ) 


Unless stated otherwise we use N\ = 30 (see Ap¬ 
pendix). This covers the range from Ai = 0.27pm to 
A 30 = 3.34pm. 

In Figure [3j we plot Q pr (A)/a e flr for different values 
of a e ff. We normalize the Q 's by a e g because accelera¬ 
tion due to radiation pressure is proportional to Q pT /a e ff 
(see Equation Q). We also show the normalized so¬ 
lar spectrum per unit ln(A). We see that the curves of 
Q P r(X)/a e ff fall off for A/a e ff > 7 r, and go to a constant 
value (proportional to 1 /a e g) for X/a e g < 1. The amount 
of acceleration due to radiation pressure per unit ln(A) 
is proportional to XB\Q(X)/a e g. For this reason, we see 
that fi will be high if the curves of Q pr (A)/a e ff are high 
near the peak of the solar spectrum. This will be true 
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Fig. 4.— /3(a e ff) f° r a solid sphere, and for the three classes of 
ballistic aggregate targets (section^. These /3 values are averaged 
over 27 grain orientations of 16 grains of each class. Error bars 
correspond to the standard deviation of the angle-averaged /3 values 
for each class 

for intermediate sized grains, as large ones will not have 
high values of Q pr (A)/a e ff due to the l/a e ff factor and 
small grains will not have high values either, because of 
the drop-off in Q for large values of X/a e g- 


5. RESULTS 

5.1. Pure Silicate Grains 

For the solar spectrum, we present calculations of 
/3(a e ff) for four different sets of geometries. In order of 
decreasing compactness, these are: a solid sphere, the 
BAM2 clusters, the BAM1 clusters, and the BA clusters, 
all with N = 32. Figure E] shows f3(a e g). The /? values 
are averaged using 27 different orientations for each of 
16 grains [J 

We see in Figure [4] that fluffiness (low values of /) 
mildly enhances /3 for a e g > 0.3/im particles, but sup¬ 
presses /3 for a e g < 0.3/im particles. We also note that 
for all models, the peak of /3 is less than unity. All the 
models converge to within a few percent in the electric 
dipole l imit (q s g —> 0) . This is in contrast to the work 


of Kohler et al. (2007), which shows grain compactness 
to have very little effect on /3 for silicate grains between 
10 -2 and 10 pg. Our res u lts are in qualitative agreement 
with Tazaki & Nomura (2015) who also find that that 
their less compact grains have lower /3 values except for 
sizes > 0.5 /zto. 

The captured Stardust particles have M ss 3 pg, cor¬ 
responding to a a ff ~ 0.6qm, assum ing a material density 


of 3.8 g/cc. Westphal et al. (2014) additionally identify 


3 We used all 16 realizations of the classes BA.32, BAM1.32, 
and BAM2.32, available at http://www.astro.princeton.edu/ 
-draine/agglom.html 


Fig. 5.— Scatter between different targets and orientations. Each 
point represents one of the 16 targets for each class at a randomly 
selected orientation. Solid curves (identical to those in Figure |4[) 
how the average for each class (BA, BAM1, BAM2) as a function 
of a e ff. In each panel, we have also plotted the value for a single 
sphere. 

impact craters on the aluminum foil which are consis¬ 
tent with much smaller particles: a e ff = 0.1 — 0.15/im. 
Both a e ff = 0.15p,m and a e g = 0.6/mr are off the peak 
in Figure |4j reducing the predicted /3 for these sizes to 
< 0.5. 

The estimated densities for the two captured Stardust 
grains were 0.7 and 0.4 g/cc, corresponding to filling fac¬ 
tors 0.1 < / < 0.2, corresponding most closely to the BA 
aggregates (see Table [I]). 

5.2. Variation of f3 with Grain Realization and 
Orientation 

The results presented in the previous section suggest 
that porous silicate grains have /3 values that peak below 
0.6. However, there is some scatter arising both from dif¬ 
ferent instances of each class of grain, and also from dif¬ 
ferent possible orientations. To investigate scatter from 
these two sources, for each grain and each size, we pick 
one random orientation, and calculate the /3 value for 
that grain, size, and angle. These points are all plotted 
in Figure 5] with the curves (imported from Figure [4| 
showing the average values for each class. This gives a 
sense of the range of likely values for /3 that we would 
see in a distribution of different grains which enter the 
solar system with different orientations. The distribu¬ 
tion shown here has a somewhat larger spread than the 
actual one, as there will be additional averaging over ori¬ 
entations inherent to whatever spin state the grains have. 
We see a scatter on the order of 10% within each class 
of grain. All of our conclusions are robust to arbitrary 
choices of grain instance and orientation angle. 

From Figure[5]we cannot tell if the majority of scatter 
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is due to variations between different grains of the same 
class, or due to different orientations. To investigate this, 
for each grain, we calculate /3 for N OT \ = 10 random ori¬ 
entations. Then for each grain class and size, we define 
two statistics 


and 


N sr 

= E 

2=1 


1 

M 


Y ' 1 {Pi,j ~ hi ) 2 

AU-1 ’ 

J= 1 


^tot 


Ngr N ori 

a EE 

M i=i j =i 


- m) 2 

N or iN gI — 1 ’ 


(9) 


( 10 ) 


where /q is the value of f3 averaged over orientations of 
a single grain, and /.i is the value of (3 averaged over all 
grains and orientations, ag is the mean standard devia¬ 
tion in f3 for a given grain sampled at different orienta¬ 
tions. ex tot is the standard deviation in /3 for a given grain 
class sampled for different instances and orientations. If 
most of the variance is due to different grain instances, 
rather than orientations, then we expect og -C cr tot . 
These quantities are tabulated for two different grain 
sizes. In the last column we compare the variance at¬ 
tributable to orientation with the total variance within 
each grain class, and come to the conclusion that the ma¬ 
jority of variation in f3 is due to grain orientation, rather 
than grain instance. 


TABLE 2 

Scatter in f 3 due to angles and instances 


Grain Class 

n tot 


(_£e_\2 
^ crtot. ' 

BA 0.15 pm 

0.0181 

0.0156 

0.743 

BAM1 0.15 pm 

0.0161 

0.0153 

0.903 

BAM2 0.15 pm 

0.0137 

0.0117 

0.729 

BA 0.6 fin i 

0.0236 

0.0201 

0.725 

BAMl 0.6 pm 

0.0187 

0.0169 

0.817 

BAM2 0.6 fim 

0.0182 

0.0170 

0.872 


5.3. Metallic Iron Inclusions 


The Stardust team reports (Westphal et al. 20141 the 
presence of “Fe-bearing phases" in two captured grains 
(“Orion” and “Hylabrook”), possibly consistent with 
metallic iron. Some interstellar grain s may contain a sub- 
stantial amount of metallic iron (see Jones|[l990 Draine 


& Hensley 2013). We calculate the effect of Fe on fj by 


replacing some of the spheres in the target structures 
with material having the density and dielectric function 
of iron, rather than astrosilicate. This was done assum¬ 
ing that the iron spheres are randomly located (that is, 
for a grain with N spheres, M of which are iron, any 


of the O configurations are equally likely to be cho¬ 
sen). Figure|6]shows a scatterplot of /3 values for different 
numbers of iron spheres and different effective radii for 
16 members of the BA class of targets. For each target, 
we picked a random orientation. We see that the addi¬ 
tion of iron substantially increases /3 for a e g < 0.3/im. 
However, to make (3 exceed one, we require > 35% of the 
solid volume be iron, and a e fi < 0.25 /im. Such grains 
have < 0.16 pg of silicate material, far below the ~ 3 — 4 
pg masses of Orion and Hylabrook. 



Fig. 6.— Curves of d vs. a e ff for targets of the BA class 
(BA.32.1-BA.32.16), with random orientations. This is done for 
targets where some of the 32 silicate spheres have been replaced 
with iron spheres of the same size. Solid lines represent averages 
of all the 16 data points. 

5.4. Carbon 


The captured grains might possibly have had carbona¬ 
ceous mantles, perhaps lost during the collection process. 
We consider three different experimentally determined 
dielect ric fu nctions for amorphous carbon from |Zubko| 


et al. (1996), and show that we cannot bring (3 above 


unity for a grain with 3 pg of silicate material by adding 
carbon to the grains. 

Figure [7] shows /3 as a function of a e ff for different val¬ 
ues of the carbon fraction. Carbon spheres were dis¬ 
tributed randomly throughout the grain as in Sectio n 


5.3 


son 


We a ssume a carbon density of 2.0 g cm 3 (Robert- 


1986). We see that f3 peaks above 1 for all three 


dielectric functions if 1/4 or more of the material by vol¬ 
ume is carbon. However, j3 > 1 occurs only for grains 
with small silicate masses. Figure [8] shows [3 for particles 
containing 3 picograms of silicate material (motivated 
by t he 3 putative interstell ar particles captured by Star¬ 
dust Westphal et al. 2014), as a function of total grain 
mass, assuming the non-silicate mass to be in the form 
of amorphous carbon. The top panel assumes the car¬ 
bon spheres to be distributed randomly throughout the 
grain, and the bottom panel assumes that the spheres 
which impact the structure later during the generation 
algorithm are carbon (leading to carbon preferentially on 
the outside of the grain). In no case do we find f3 > 1 for 
grains with 3 pg of silicate material. 


5.5. Number of Spheres 

We also consider changing the number of spheres in our 
model dust grain. We consider model grains generated 
by the same algorithm, but with N = 256 (the largest 
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Fig. 8.— /3 as a function of total mass for grains which contain 3 
picograms of silicate material. The upper panel shows grains with 
N = 32 spheres, with carbonaceous spheres distributed randomly. 
The lower panel shows grains with N = 32 spheres, where the 
silicate spheres arrive first during the formation process, and the 
carbon spheres arrive later, forming a fluffy “mantle”. Error bars 
represent the 1 — cr scatter from grain to grain. 

N for which the calculations remain tractable) instead 
of 32. Figure [9] is analogous to Figure [5] for each mem¬ 
ber of each target class, we pick a random orientation, 
and calculate /3. We see that there can be substantial 
differences between N = 32 and N = 256 (the filling 
factor drops by ~ 20% as N increases from 32 to 256), 
particularly in the value of a e g where /3 peaks, but our 
principal conclusion — that /3 peaks substantially below 
unity for silicate clusters — remains robust regardless of 
the number of spheres used. 


Fig. 7.— Each panel of this figure is analogous to Figure [6] for 
a different carbon dielectric function. 
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a eff 


Fig. 9. — Variation of 0 depending on number of spheres in the 
target cluster. Each point represents one of the 16 targets for 
each class at a randomly selected (from an isotropic distribution) 
angle. Each panel represents a different grain model, as indicated. 
Instances with N = 32 are represented with filled circles. Instances 
with N = 256, are represented by crosses. Solid lines represent 
average values. These are taken from Figure[4]for the case with 32 
spheres, and calculated as the mean of the data points for the 256 
spheres case. 


6. TRANSVERSE FORCES 

When modeling the dynamics of dust grains entering 
the solar system, it may not be appropriat e to ig nore the 
trans verse force from radiation pressure (Kimura et al. 
2002). It is difficult to correctly calculate the transverse 
force averaged over the rotation of the grain, as the spin 
state depends on the initial spin state far from the Sun, 
torques from solar radiation pressure, and internal relax¬ 
ation processes. Here, we present an estimate of the im¬ 
pulse delivered to the grain due to transverse radiation 
forces. For each of the target geometries, we used the 
discrete dipo l e app roximation code DDSCAT (Draine & 


Flatau 1994 2013). Because DDSCAT is slower than 
the MSTM coae^ we used N\ = 10 instead of 30 to speed 
up calculations. Based on the results in the appendix, 
we expect this to make less than a 2 % difference to the 
results. 

We calculated values of the “transverse p” or p t . This 
is the ratio of transverse radiation force to gravitational 
force, pt will vanish if we average over all grain orienta¬ 
tions, so we calculated it for grains spinning about their 
principal axis of largest moment of inertia, for different 
values of the angle 9 a \ between the spin axis and the di¬ 
rection to the Sun. This was calculated for a e g = 0.15 
and 0.6 0 m grains of the BAM1 geometry. Results are 
shown in Table [3] Because we average over rotations 


about the spin axis, Pt(9 a i = 0) = 0 by symmetry. 


TABLE 3 

Transverse 0 statistics for BAM1, N = 32 (realizations 
BAMI.32.1 - BAMI.32.16) 


a eff j 

8 a 1 

(ft) 

max(/3t) 

0.15 

0° 

0 

0 

0.15 

CO 

O 

o 

0.067 

0.098 

0.15 

o 

O 

CD 

0.087 

0.129 

0.15 

CD 

O 

o 

0.0051 

0.0097 

0.6 

0° 

0 

0 

0.6 

CO 

o 

o 

0.011 

0.026 

0.6 

o 

O 

CD 

0.021 

0.047 

0.6 

90° 

0.0077 

0.020 


We can estimate the transverse velocity arising from 
radiation pressure by considering an instructive case with 
P = 1, and Pt <C 1. This case has straight-line constant- 
velocity particle trajectories except for a small pertur¬ 
bation due to pt. Let us assume that the transverse 
force acts in the same direction over the trajectory. This 
would be the case if both the transverse force and the 
grain’s spin axis were perpendicular to the orbital plane, 
but provides a reasonable estimate and an upper bound 
for more complicated geometries. In this case the total 
transverse velocity acquired by the grain on its journey 
from infinity to the point ( R , 9) (where R is the grain- 
Sun distance, and 9 is the angle between the interstellar 
wind direction and the grain-Sun vector) is: 

v= r PtGM Q dx PtGMp 9 
Vt Jrcos 6 vpx 2 + (Rsin9) 2 ) Rvi sin6>’ 


where Vi is the velocity of the Sun with respect to the 
local ISM. For numbers (R = 2AU, = 26 km s -1 and 


9_== 60°) appropriate for the Stardust mission (Sterken 
et al. 2014), this gives v t = 21/3 t km s -1 . Using typical 
numbers from Table f3j we see that we expect changes in 
transverse velocity of~ 1 km s -1 for the small grains that 
left impact craters in the Stardust mission, and ~ 0.2 km 
s _1 for the large captured grains. Including a true value 
of p less than one should slightly lower the transverse 
velocity, as the transverse force has less time to act due 
to the increased radial velocity. 


7. TORQUES 

We can also compute radiative torques on dust grains 
using the DDSCAT code. We computed the component 
of these torques along the spin axis for the ‘BAM1’ class 
of grains, with effective radii of 0.15 and 0.6 /irn. As 
in the previous section we calculated torques for grains 
spinning about their principal axis of largest moment 
of inertia, for different values of the angle 9 a i between 
the spin axis and the grain-Sun vector. We evaluate the 
torque using DDSCAT. The torque efficiency Qr(A) is 


Qr = 


2ttF 


^H e ff Hrad^ 


( 12 ) 


(Draine & Weingartner 1996), where T is the torque, and 
■u rac j the energy density of radiation with wavelength A. 











































Given a spectrum of flux F( A), the torque is 
f = Tral s u ia . d ^(Q r ), 


(13) 


where M ra d is the energy density of the incident radiation, 


and 


. _ / 0 °° XF(X)Qrd\ 

Wr) f 0 °°\F(\)d\ 

£aw 

1 /o°° F(X)d\ 


(14) 


(15) 


Assuming the trajectory to be undeflected by the combi¬ 
nation of radiation pressure forces and solar gravity (i.e. 
/3 = 1, /3 t = 0 ), we can integrate Y(r)dt = Y{r)dr/v r {r) 
along the grain trajectory (analogous to the calculation 
in Section pi , assuming zero spin angular momentum at 
infinity, to nnd the spin angular momentum of the grain: 


L(R) 


(Qr)a 2 eB L Q (X) 0 

8 TTViRc sin 0 


(16) 


This calculation neglects rotational damping due to 
emission of infrared photons. The damping torque on 
a spinning grain from radiation of photons with wave- 
len gths long compared with t he size of the grain is given 
by ( Ali-Haimoud et ah1|2009 ) 


dL _ -ui r°° P„dv 

dt 2ir 2 J 0 v 2 


(17) 


where P v is the radiated power per unit frequency. Let 
the grain have temperature T gr = yTee, where Tbb = 
To (0.5 Rq/R) is the “blackbody” temperature. Submi¬ 
cron grains are poor thermal radiators, and will likely 
have 7 > 1. For an infrared opacity oc v s , the damping 
torque is 


dL 

dt 


-2 

7T7 2 


T (2 + s)C(2 + s) h 2 
T(4 + s)C(4 + s) k 2 


(T ^BB( ( 3abs)ao ff a;, 


(18) 


where (Q a bs) is the absorption efficiency averaged over 
the solar spectrum, and £( x ) is the Riemann £ function. 
For s = 2, this becomes 


— —0.034 ^ 2^2 ^BB^absVeffW. 

Since 

the damping timescale t = L/L is 

T = _ 87rpq^ ff fc 2 7 2 _ 

15 x 0.034/i 2 crT| B (<5 abs )/ 2 / 3 

27yr 2 ( « eff \ 3 / R \ 2 /O^A 273 
(Q abs ) 7 \0.2/rm/ \2AU/ \f) 


(19) 

( 20 ) 


( 21 ) 


Given the expectation that 7 > 1, this is much longer 
than the timescale on which grains gain their spin (which 
is just the dynamical time R/v ~ 0.37i?/(2AU) years, 
since torque is proportional to the gravitational acceler¬ 
ation of the grain, so L should vary on the same timescale 


as v) for 0.15-0.6 [i m grains at 2 AU, so we can ignore 
spin damping. 

Using Equation (fl6j) and letting 0/ sin0 ss 1, Table [4] 
shows the mean and maximum surface velocity u sur f « 
La e ff/^ 173 // for a e ff = 0.15/im and 0.6 /jm grains as a 
function of 9 a 1 , the angle between the rotation axis and 
the direction to the star. These are calculated from an 
ensemble of 16 grains of each size. 


TABLE 4 

Surface velocity for interstellar BAM1 grains at R « 2 

AU 


®eff 5 f-LTCl 

dal, degrees 

mean ^ SU rf? ms 1 

max t?surf ? m, m s 1 

0.15 

0 

458 

945 

0.15 

30 

368 

727 

0.15 

60 

205 

495 

0.15 

90 

74 

226 

0.6 

0 

471 

1473 

0.6 

30 

230 

560 

0.6 

60 

72 

197 

0.6 

90 

115 

241 


Having estimated surface velocities, we then ask 
whether this will lead to centrifugal disruption. We can 
estimate the critical velocity by equating the centrifu¬ 
gal force pulling two hemispheres apart with the cross- 
sectional area times the yield stress S ma *. This calcu¬ 
lation shows that centrifugal disruption will occur for 
v > u cr it = 2 yjS max /p. For ideal materials w ith no frac¬ 
tures, S rn „ ~ 10 11 dyn cm -2 (MacMillan 19721, and 


u cr ;t is of order the sound speed in the material. 

However, because of defects, real materials have tensile 
strengths well below the ideal value. If we use a density 
of 2.4 g cm -3 , and S ma , x = 5.0-10 7 d yn cm~ 2 , app ropriate 
for construction grade concrete (Anoglu et al.j 20061 , we 
obtain v cr it ~ 90 m s” 1 . This is smaller than the esti¬ 
mates for v SUI f in Table[4| implying that some grains may 
be centrifugally disrupted by the solar radiation torques. 

The spin-up of submicron grains by solar radia¬ 
tion is analogous to the YORP (Yarkovsky- O ’Keefe- 
Radzievskii-Paddack) effect (see Vokrouhlicky et al. 
(20151 for a review), but where the irregularities are 


on scales comparable to the wavelength. Dust grains 
have far higher critical angular velocities for disruption 
than rubble piles because the former are held together 
by chemical bonds rather than gravity. 


8. PARTICLES CAPTURED BY STARDUST 

The Stardust mission captured two particles (“Orion” 
and “Hylabrook”) that are claimed to be consistent with 
interstellar dust particles entering the solar system. The 
condition of these particle s is consistent with impact ve¬ 
locities <C 10 km s _1 (Westphal et al. 2014). The particle 
masses correspond to a e ff ~ 0 . 6 pm, assuming a material 
density ~3.8 g cm -3 . 

Impact speeds <10 kms -1 require /3 > 1, so that so¬ 
lar radiatio n pressure can deceler ate incoming grains (see 
Figure 8 of Sterken et al. p014|). Using silicate grains, 
supplemented with carbon and iron, we are unable to 
find a composition such that /? > 0.9 for a gra in with 
3 picograms of silicate material (see Figures [ 6 ] and [ 8 ]). 
Therefore, if the actual impact speeds are indeed < lCkni 
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s -1 , it seems unlikely that particles Orion and Hylabrook 
originate in the interstellar stream approaching the solar 
system at ~ 26 km s -1 ; some other origin, presumably 
in the solar system, appears to be required. While we 
have not considered every possible set of grain proper¬ 
ties, we have covered much ground, both in composition 
and structure, without managing to produce a grain with 
high enough (3. 

9. SUMMARY 

1 . Accurate calculations of models for “fluffy” grains 
with the dielectric functions and density appropri¬ 
ate for astronomical silicates show that /3 < 1 for 
all sizes. High grain porosity does not enable sub¬ 
micron sized grains to have substantially higher 
ratios of radiation pressure to mass. In fact, for 
a e ff = 0.1 — 0.15/rm grains (such as those responsi¬ 
ble for the Stardust A1 foil impact craters), higher 
porosity tends to substantially reduce (3 at constant 

®eff • 

2. One way to potentially increase /3 is to include 
metallic iron. However, we find that [3 peaks below 
one unless about half the mass of the grain is iron. 
For the grains with the masses of those captured 
intact in aerogel by the Stardust mission, /3 would 
be < 0.6 even if the grain were entirely iron. 
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3. We also consider mixed silicate-amorphous carbon 
structures. For grains with ^-3pg silicate mass (as 
for the Orion and Hylabrook grains captured by 
Stardust) we find /3 < 0.8 (see Fig [ 8 ]). For par¬ 
ticles with mass and composition resembling the 
captured Stardust grains, we are unable to find a 
single example with large enough /3 for the particle 
to be of interstellar origin, even if they are allowed 
to have amorphous carbon mantles that were lost 
during capture. 

4. Transverse forces from radiation pressure on incom¬ 
ing interstellar grains lead to velocity changes of at 
most ~2 km s _1 at 2 AU. 

5. Radiative torques due to sunlight can drive irregu¬ 
lar sub-micron grains entering the solar system to 
spin with equatorial velocities of a few hundred me¬ 
ters per second. Depending on the tensile strength 
of the grains, this could lead to rotational break-up. 
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APPENDIX 

We investigate the sensitivity of the calculated (3 to the choice of N\, the number of wavelengths sampled. Figure 
[T0| shows the magnitude of the fractional change in (3, relative to the result for N\ = 30, for different grains as labelled 
on the plot. This is averaged over one orientation of 16 different grains in each class. We see very good convergence 
at the few percent level for all N\ > 5. 
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Fig. 10.— Magnitude of fractional change in ft when using number of wavelengths N\ < 30. Errors are under 5% relative to N\ = 30 for 
values of N\ as low as 5. 











